UNCLASSIFIED 


_ AD  NUMBER _ 

AD848843 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  DoD  only; 
Administrative/Operational  Use;  SEP  1968.  Other 
requests  shall  be  referred  to  Space  and  Missile 
Systems  Organization,  Los  Angeles  AFB,  CA. 


_ AUTHORITY 

SAMSO  per  DTIC  form  55 


THIS  PAGE  IS  UNCLASSIFIED 


SAMSQ-TR-68-439 


09601-60 J6-T0-00 


TRANSFORMATION  THEORY  FOR  THE  COMPRESSIBLE  TURBULENT 
umiNllARY  LAYER  WITH  ARBITRARY  PRESSURE  GRADIENT 


by 

John  E.  Lewis,  Toshi  Kubota,  and  Wilmot  H.  Webb 


September  1968 


Advanced  Research  Projects  Agency 
ARPA  Order  No.  888 
Contract  No.  F04701-68-C-0G41 


Prepared  for 

Space  and  Missile  Systems  Organization 
Air  Force  Systems  Command 
Norton  Air  Force  Base,  California  92409 


rj  Q 

siriPfaCTonra 

MAR  14 1969 

[ISSEirU'El 


Each  transmittal  of  this  document  outside  the  Department  of 

Defense  must  haYe  prior  approval  of.  Headquarters, ^^and 

Missile  Systems  Organization  ( 

Also,  this  document  is  subject  to  special 

expor ^controls  and  each  transmittal  to  foreign  governments 
or  foreign  nationals  may  be  made  only  with  prlor^approv^^^ 
of  Space  and  Missile  Systems  Organization  {jSmmammaMOkim 

samso  i  msd 

AFS  CALIF  90045 


TRltf 

SYSTEMS  GROUP 


f-  ,  \ 

A.;i>  ^ 

SAMS0^TR-68-439 

^ _ ^  I 


096O1-6036-T0-O0 


\ 

TRANSFORMATION  theory  for  the  compressible  turbulent 
BOUNDARY  LAYER  WITH  ARBITRARY  PRESSURE  CRADTENT. 


0  !  *  )  x  c  c  \  k-  T  Li  t  .  y 

1 -  r 


by 


^  John  E./  Lewis  ,  ToshiyKubota ^  Wilmot  H./Webb 


jV/  SeptWMtor  £968 


Advanced  Research  Project 
ARPA  Order  No.  888 
Contract  No.l  FJ047^01-68- 


Prepared  for 

Space  and  Missile  Systems  Organization 
Air  Force  Systems  Command 
Norton  Air  Force  Base,  California  92409 


s  Agency  / 


Each  transmittal  of  this  document  outside  the  Department  o^ 
Defense  must  have  prior  approval  of  Headquarters,  Space  an 

Missile  Systems  Organization  ( 

Also,  this  document  is  subject  to  special 

export  controls  and  each  transmittal  to  foreign  governments 
or  foreign  nationals  may  be  made  only  with  prior  approval 
of  Space  and  Missile  Systems  Organization 

ill  *****  |  rfTJrff,l‘"  ,  D  o:  >OU 

L.A.  AFS  CAL! I  90041 


Prepared  for 

SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
AIR  FORCE  SYSTEMS  COMMAND 
NORTON  AIR  FORCE  BASE 
CALIFORNIA  92409 


Approved  by 


. 

L.  A.  Hronas 


Manager 

Fluid  Physics  Department 


TRW  Systems  Group 
Redondo  Beach,  California 


ii 


FOREWORD 


This  research  was  supported  by  the  Advanced  Research  Projects  Agency 
of  the  Department  of  Defense  and  was  monitored  by  the  Space  and  Missile 
Systems  Organization  of  the  Air  Force  Systems  Command,  Norton  Air  Force 
Base,  San  Bernardino,  California.  The  work  reported  here  is  submitted  as 
a  technical  report  on  Task  7.4,  Chemically  Reacting  Turbulent  Boundarv  Layer, 
of  Penetration  Systems  Studies,  Contract  F04701-68-C-0041. 

Each  transmittal  of  this  document  outside  the  Department  of  Defense 
must  have  prior  approval  of  Headquarters,  Space  and  Missile  Systems  Organ¬ 
ization  (SMYSE),  Norton  Air  Force  Base,  California  92409.  Also,  this 
document  is  subject  to  special  export  controls  and  each  transmit  taj.  to 
foreign  governments  or  foreign  nationals  may  be  made  only  with  prior  ap¬ 
proval  of  Space  and  Missile  Systems  Organization  (SMYSE),  Norton  Air  Force 
Base,  California. 

The  above  statement  is  made  for  the  reason  that  information  in  this 
report  is  embargoed  under  the  Department  of  State  International  Traffic 
in  Arms  Regulations.  This  report  may  be  released  to  foreign  governments 
by  departments  or  agencies  of  the  U.  S.  government  subject  to  approval  of 
Space  and  Missile  Systems  Organization  (SMYSE),  Norton  Air  Force  Base, 
California,  or  higher  authority  within  the  Department  of  the  Air  Force. 
Private  individuals  or  firms  require  a  Department  of  State  export  license. 

Publication  of  this  report  does  not  constitute  Air  Force  approval  of 
the  report's  findings  or  conclusions.  It  is  published  only  for  the  exchange 
and  stimulation  of  ideas. 

Major  Walter  D.  McComb/SMYSE 
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ABSTRACT 


We  have  rex iewed  Coles'  analytical  development  of  the  transformation 
and  its  application  to  the  constant  pressure  turbulent  boundary  layer. 

The  limitations  of  the  theory  noted  by  Baronti  and  Libby  are  clarified, 
and  an  attempt  is  made  to  explain  the  remaining  discrepancies.  The 
equations  governing  the  scaling  functions  (which  link  the  high  speed 
boundary  layer  to  an  equivalent  low  speed  flow)  are  derived  in  terms  of 
the  integral  properties  of  the  low  speed  flow  for  an  arbitrary  pressure 
distribution.  The  modeling  is  completed  by  coupling  the  pressure  gradients 
of  the  two  flows  and  predicting  the  behavior  of  the  integral  properties 
of  the  low  speed  flow.  The  resulting  formulation  predicts  the  integral 
properties  of  the  high  speed  turbulent  boundary  layer  as  well  as  the 
required  modeling  to  an  equivalent  low  speed  flow.  Finally,  the  theory 
is  applied  to  a  variety  of  situations. 
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1.  INTRODUCTION 


The  dynamic  equations  governing  the  boundary  layer  are  identical 
in  form  for  laminar  and  turbulent  flow  with  the  exception  of  the  specific 
form  for  the  shear  stress.  While  the  simple  Newtonian  stress  relation¬ 
ship  is  valid  everywhere  for  laminar  flow,  it  is  known  to  be  valid  only 
very  near  the  wall  in  a  turbulent  boundary  layer.  The  prediction  of  the 
turbulent  boundary  layer  is,  of  course,  much  more  difficult  because  of 
the  lack  of  understanding  of  the  turbulent  momentum  transport.  Most 
analytical  models  rely  on  a  simple  gradient  diffusion  process  (eddy 
viscosity); and  while  such  a  method  is  perhaps  adequate  for  low  speed 
flows,  the  effect  of  compressibility  on  such  a  model  is  unknown.  This 
difficulty,  in  addition  to  the  existence  of  a  vast  amount  of  empirical 
knowledge  about  the  low  speed  turbulent  boundary  layer,  makes  the  prospect 
of  using  a  mathematical  transformation  to  account  for  the  effects  of 
compressibility  particularly  inviting. 

The  objective  of  the  present  study  is  to  develop  a  transformation 
theory  which  will  model  a  given  supersonic  turbulent  boundary  layer  flow 
of  arbitrary  pressure  gradient  [u  (R  )  -  specified]  into  an  equivalent 
low  speed  boundary  layer  flow  [u  (R— )  -  predicted].  It  will  be  found 
that  in  developing  this  modeling,  the  integral  properties  of  both  the 
high  speed  and  the  equivalent  low  speed  flow  will  be  predicted  as  well. 

We  begin  by  reviewing  Coles’ ^  analytical  development  of  the 
transformation  and  its  application  to  the  constant  pressure  turbulent 
boundary  layer.  The  inaccuracies  of  the  theory  noted  by  Baronti  and 
Libby v  '  are  clarified,  and  an  attempt  is  made  to  explain  the  remaining 
discrepancies . 

The  equations  governing  the  scaling  functions  (which  link  the  high 
speed  boundary  layer  to  an  equivalent  low  speed  flow)  are  derived  in 
terms  of  the  integral  properties  of  the  low  speed  flow  for  an  arbitratry 
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pressure  distribution.  The  modeling  is  completed  by  relating  the  pressure 
gradients  of  the  two  flows  and  predicting  the  behavior  of  the  integral 
properties  of  the  low  speed  flow.  The  resulting  formulation  predicts  the 
integral  properties  of  the  high  speed  turbulent  boundary  layer  as  well 
as  the  required  modeling  to  an  equivalent  low  speed  flow. 

Finally,  the  theory  is  applied  to  a  variety  of  situations. 
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2.  REVIEW  OF  COLES'  TRANSFORMATION  WITH  CONSTANT  PRESSURE  APPLICATIONS 


The  use  of  a  mathematical  transformation  to  account  for  the  effect  of 
compressibility  has  been  extensively  applied  to  laminar  boundary  layers  and 
wakes.  Howarth, ^  Stewartson, ^  and  many  others  have  used  this  technique. 
Coles  has  formulated  perhaps  the  most  general  transformation  yet  considered 
for  the  dynamic  boundary  layer  equations.  These  equations  for  the  steady, 
two-dimensional  mean  flow  of  a  compressible  fluid  are 

h (pu)  +  h (cv>  * 0  (1) 


pu 


_  +  ii 

dx  dy 


(2) 


The  transformation,  simply  stated,  specifies  the  correspondence  between 
a  flow  governed  by  these  equations  and  an  incompressible  flow  for  which 


9u  t  3v 
dx  dy 


0 


(3) 


- 8u  . - 9u 

P  U  “  +  p  v  — 

dx  dy 


.  ^£  +  2i 

dx  dy 


(4) 


and  p  is  constant. 

The  transformation  set  forth  by  Coles  is  represented  by  thr-  e 
initially  unspecified  scaling  functions: 


The  term  transformation  will  be  used  here  in  the  same  sense  used  by  Coles, 
i.e. ,  as  a  true  mapping  of  one  flow  field  into  another  rather  than  simply 
a  mathematical  manipulation. 
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dx 

dx 


(6) 


-  e 


p_  jx 

p  9y 


n 


(7) 


In  order  to  insure  the  boundedness  of  the  pressure  gradient  throughout 
each  flow  field,  Coles  imposes  the  condition  that  these  functions  be 
independent  of  y  ,  i.e.,  a  -  a(x),  n  ■  n(x)  and  £  -  £(x)  .  Before 
pursuing  this  development,  it  is  appropriate  to  comment  on  the  ability 
of  the  transformation  to  predict  the  behavior  of  high  Mach  number  boundary 
layers  at  constant  pressure. 

Baronti  and  Libby  in  attempting  to  evaluate  the  transformation  at 
nominally  constant  pressure  noted  a  systematic  deviation  with  increasing 
Mach  number  in  the  "wake"  region  of  the  boundary  layer  but  concluded 
that  the  "law  of  the  wall  region"  was  well  predicted. 

They  arrive  at  this  conclusion  by  choosing  a  value  of  Cf  which 
allows  the  best  agreement  of  the  respective  low  speed  and  high  speed 
velocity  profiles  in  the  inner  logarithmic  region  of  the  boundary  layer. 
However,  this  procedure  precludes  a  test  of  the  transformation  in  the  wall 
region  of  the  boundary  layer  since  it  forces  a  match  there;  in  addition, 
it  does  not,  in  general,  satisfy  an  essential  requirement  of  the  trans¬ 
formation  which  relates  one  station  to  another  in  the  high  speed  and  low 
speed  flows,  i.e.,  Coles'  "law  of  corresponding  stations",  which  for 
y  r~  T  is  simply  CfRQ  -  CfR-^  .  An  additional  but  less  fundamental  point 
is  that  the  choice  by  Baronti  and  Libby  of  6  ,  the  boundary  layer  thickness, 
for  normalization  gives  rise  to  uncertainties  which  can  be  large 
particularly  for  turbulent  boundary  layers  where  the  edge  is  not  well 
defined. 

Using  the  law  of  corresponding  stations,  C^R^  ■  C^R—  ,  and  comparing 
the  velocity  profiles  in  the  coordinates 
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u 

u 


u 

u 


vs. 


1 

"o 


p_ 

p. 


dy 


0 


we  eliminate  both  of  the  above  objections  to  the  Baronti-Libby  comparison 

and  the  need  to  specify  any  of  the  scaling  functions,  hence,  the  need  to 

resort  to  any  ad  hoc  theory  such  as  the  "sublayer"  or  "substructure"  hypothesis. 

The  results  of  applying  the  transformation  in  this  way  to  the 

experimental  data  of  Coles  ^  and  Korkegi^  (adiabatic,  flat  plate)  is 

(7) 

shown  in  Figure  1  contrasted  with  the  low  speed  data  of  Wieghardt 

(1)* 

and  the  analytical  form  suggested  by  Coles.  As  can  be  seen,  there  is 

indeed  a  deviation  from  the  low  speed  profiles  which  increases  as  Mach 
number  increases.  However,  the  deviation  is  seen  to  be  more  significant 
in  the  inner  region  than  the  outer  region  of  the  boundary  layer  which  is 
contrary  to  the  conclusions  of  Baronti  and  Libby.  In  addition,  the  thickness 
of  the  boundary  layer  is  increased.  The  actual  magnitude  of  the  velocity 
in  the  wake  region  is  within  2-3%  of  the  predicted  value  whereas  the  deviation 
at  y/0  ■  1  is  seen  to  be  as  large  as  10%. 

Apparently,  then,  the  transformation  fails  in  the  inner  "law  of  the 

wall"  region  of  the  boundary  layer  to  an  extent  which  increases  systematically 

with  Mach  number.  Since  this  region  typically  represents  less  than  20% 

of  the  boundary  layer  thickness  and  it  can  reasonably  be  argued  that  the 

resulting  error  in  overall  (e.g.,  integral)  properties  of  the  layer  may 

not  be  excessive,  nevertheless  the  fundamental  idea  of  Coles  that  a  simple 

geometric  mapping  of  a  low  speed  to  a  high  speed  flow  can  account  for 

effects  of  compressibility  does  not  appear  to  be  valid  in  detail.  One 

phenomenon  which  is  not  expected  to  be  accounted  for  by  simple  transformation 

theory  is  the  increasing  magnitude  of  sound  pressure  level  [the  irrotational 

(8)  (9) 

fluctuation  mode  in  the  parlance  of  Kovasznayv  1  and  Morkovinv  ']  with 


The  analytical  form  has  the  ambiguity  previously  mentioned  regarding  6  ; 
however,  here  it  is  confined  to  an  analytical  form  and  affects  only  one 
of  the  two  flows. 
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Mach  number.  Whether  this  phenomenon  is  important  to  the  momentum  exchange 
is  not  clear  [see  Schubert  and  Corcos^10^],  but  it  is  difficult  to  see  how 
the  complicated  sound  radiation  field  generated  at  high  Mach  number  in  the 
boundary  layer  could  be  accounted  for  merely  by  a  mathematical  transformation. 
Another  phenomenon  introduced  by  high  Mach  number  is  the  variation  of 
molecular  transport  coefficients  (viscosity,  conductivity)  due  to 
temperature  variations  in  the  boundary  layer.  An  important  process  affecting 
turbulent  momentum  exchange  in  the  region  adjacent  to  the  laminar  sublayer 
is  the  "damping"  of  velocity  fluctuations  by  viscous  forces  there.  For 
an  adiabatic  supersonic  boundary  layer,  the  temperature  near  the  wall  is 
significantly  higher  than  ambient  and  an  increased  dissipation  of  turbulent 
velocity  fluctuations  due  to  the  increased  viscosity  is  therefore  expected 
in  this  region.  This  increase  in  viscous  dissipation  is  believed  to  be  a 
plausible  explanation  for  the  deviation  of  the  transformed  profiles  in  the 
law  of  the  wall  region.  To  give  some  quantitative  support  to  this  conjecture, 
we  have  made  a  rough  estimate  of  this  effect  in  the  following  way.  Webb 
has  arrived  at  an  empirical  estimate  of  the  effect  of  dissipation  based  on 
a  simple  gradient  diffusion  model  which  seems  to  give  an  accurate  representa¬ 
tion  of  the  turbulent  shear  stress  near  the  wall  in  confined  low  speed 
turbulent  flows.  This  estimate  is  obtained  by  applying  the  known  solution 
for  the  decay  of  velocity  in  a  finite  viscous  vortex  to  the  estimate  of  the 
damping  of  the  turbulent  velocity  fluctuations  which  are  convected  during 
the  mixing  process.  The  result  is  that  the  local  shear  stress  is  reduced 
by  a  factor  dependent  on  the  local  ratio  of  turbulent  to  laminar  shear 
stress,  T'  ■  t  /t  .  Figure  2  shows  the  change  in  velocity  profiles  for 
low  speed  pipe  flow  which  results  from  arbitrarily  increasing  the  value 
of  the  dissipation  factor  by  various  amounts.  Qualitatively,  the  change 
appears  similar  to  that  due  to  increasing  Mach  number  in  the  compressible 
adiabatic  boundary  layer. 

A  typical  value  of  the  ratio  Tf  ,  evaluated  at  y  -  6  is  estimated 
for  a  compressible  boundary  layer  as  follows: 
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T 1 


_p_82(3u/jy)2^iL  /e)__e_-Jj7  yg-  \  /M 

U  3u/3y  vw  T  vu  V  pw  o  Cf  ,  L) 


1/2 


Taking  P  ~  T  and  p ~  T  and  recalling  that  CfRQ  «  C  R-  ,  then  the 
ratio  of  T  for  low  speed  flow  to  its  valu  for  compressible  flow  is 


To  see  whether  the  parameter  indeed  collapses  the  velocity  profiles 
for  various  Mach  numbers  and  Reynolds  numbers,  data  taken  by  Coles ^  and 
by  Lobb,  Winkler  and  Perschv  over  a  range  of  M^  up  to  8.2  (including 
heat  transfer)  for  J  C?  6  are  shown  in  Figure  3.  The  remarkable  agreement 
of  these  profiles  for  a  fixed  value  of  “3”  strongly  supports  the  offered 
explanation  for  the  failure  of  the  transformation  theory  near  the  wall. 

As  an  alternate  display  of  the  correlation,  the  maximum  deviation  of  the 
experimental  profiles  from  the  transformation  theory  is  plotted  directly 
against  the  parameter  "J*  in  Figure  4.  Shown  for  comparison  is  the 
calculated  deviation  for  pipe  flow  referred  to  earlier. 

On  the  basis  of  these  correlations,  we  tentatively  conclude  that 
viscous  dissipation  effects  are  responsible  for  the  deviation  of  the 
transformation  theory  near  the  wall  and,  at  least  in  principle,  the 
transformation  theory  is  invalid.  However,  for  not  excessively  high  Mach 
numbers,  the  deviation  is  small,  it  extends  over  only  a  small  fraction  of 
the  boundary  layer,  and  it  does  not  preclude  the  application  of  the  theory 
to  obtain  overall  results. 
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3.  GENERALIZATION  TO  VARIABLE  PRESSURE  ADIABATIC  WALL 

3.1  SCALING  FUNCTIONS 

For  the  constant  pressure  turbulent  boundary  layer,  Coles  derives  the 
following  relationships  for  the  scaling  functions  o,  n,  S  : 


£ 

n 


U00 

—  (constant) 
u 

00 


(8) 


L  B  d(q6) 

n - d0 

P  b 


(9) 

(10) 

(y/T  ■  constant) 


He  arrives  at  these  equations  by  invoking  the  condition  of  constant  pressure, 
Newtonian  friction  at  the  wall,  and  an  assumption  known  as  the  "substructure 
hypothesis"  (which  is  similar  to  the  constancy  of  the  sublayer  Reynolds 
number).  The  analysis  is  well  detailed  in  Coles  report  and  will  not  be 
repeated  here. 


If  we  take  the  substructure  hypothesis  to  be  independent  of  pressure 
gradient  and  heat  transfer  [see  Coles ^],  we  need  only  generalize  equations 
8  and  9.  This  assumption  has  been  made  in  the  following  developments. 

While  plausible,  it  is  experimentally  unverified. 


The  pressure  gradients  of  the  two  flows  are  related 


(1) 


by 


d£ 

dx 


p  /  1  dp  +  Ue  d(n/q) 

„  2  Ip  dx  n To  dx 
5  u  \  e 


(ID 


This  equation  is  one  of  the  required  relationships;  however,  it  is  in  a 
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form  which  is  not  useful  for  our  purposes  in  that  it  involves  both 
dp/dx  and  dp/dx  ,  one  of  which  is  unknown  a  priori.  However,  it  can  be 
shown  that  by  noting  the  behavior  of  the  velocity  in  the  vicinity  of  the 
wall  [see  Lewis^13^]  the  pressure  gradients  are  linked  by  the  following 
relationship: 


_e_  _6_  dp 

p  t  dx 
w  w 


6  dp 

t  dx 
w 


(12) 


Combining  equations  11  and  12  after  some  manipulation,  we  find  that 


1 

n/o 


dn/o 

dx 


Pw  pw  pw  anl  1  dp 
Pe  py  I  P  Ug  dx 


(13) 


where,  for  reasons  which  will  become  obvious  later,  o,  ri ,  £  are  taken  as 
functions  of  x  rather  than  x  .  Equation  13  replaces  the  constant  pressure 
condition  n/o  -  (u^/  uj  (constant). 

One  more  equation  is  required  and  for  this  we  utilize  the  relationship 

(1) 

between  the  shear  stress  terms  of  the  momentum  equation, 


ii  p  g2  fl(H  11 H  (_L  _  I\  +  _i_  dn/p  L2  _  U2Y1 

-  “  2  [p  \3y  a  9y  d xj  dx  \p  P J  n/o  dx  \  e  /J 


3y  5  n 


(14) 


Integrating  across  the  boundary  layer,  we  find  that 


in  7  .  t  +  • 

2  w  w  a  dx 
a 


i  j 

*  37  dy 


+  42-  • 

dx 


(1  -  p/pe)  dy  -  '  P-“!  • 


e  e 


(1  -  u2/u2)  -p-  dy 
Me 
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Noting  that 


(1  _  u2/u2)  -2-  dy 
e  o 


5 


(o  +  0) 


i 

and  evaluating  the  integral,  j  (1  -  c/p  )  dy  by  approvim, 

(  j  ^ 

temperature  field  by  the  equation  (T  /T  )  -  1  [(T0  /: 

uirive  at  the  required  relationship 


c  u 

w  w  o  n 

~T  “  C 


Cf/2  -  07 r)  %  +  -A_  (1  +  H) 

dx  o  u*"  dx 
e 


Cf/2 


p  u  dx 
e 


We  have  tor  convenience  redefined  the  scaling  functions  in 
manner : 


l’ 


U  /  L 

>  CO  00 


The  resulting  equations  become 


1 

r'/o' 


dr  '  /  o 

dR- 

X 


(1 


+ 


u  dR— 
e  x 
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iting  the 

w)  -  U  ^ 


(1 -) 


the  following 


Cl™) 


(16,i 


) 


(19) 


Hence,  given  a  low  speed  flow  which  is  completely  specified  (at 
least  insofar  as  its  integral  properties  are  concerned;  these  three 
relationships  can  be  used  to  construct  an  equi/alent  supersonic  flow 
for  a  given  . 

The  objective  of  the  study  is,  of  course,  to  take  a  specified 
velocity  distribution  in  the  supersonic  flow,  ue(Rx)  »  and  predict  the 
necessary  u  (R— )  to  produce  its  low  speed  equivalent. 

6  X 

Returning  to  equation  (12) ,  it  can  be  shown  that 


u  dR— 
e  x 


w 


(1  +  me) 


oTnT 


u  dR 
e  x 


(20) 


Hence,  given  u  (R  )  we  have  an  equation  to  make  such  a  prediction 
provided  that  o'  and  n'  are  known.  These  functions  can  be  predicted 
by  simply  calculating  the  integral  properties  of  the  low  speed  flow 
simultaneously.  Coupling  such  a  calculation  with  the  previously  derived 
transformation  laws,  we  have  then  arrived  at  a  self-contained  set  of 
equations  which  predicts  the  modeling  we  seek.  The  integral  properties  of 
the  supersonic  turbulent  boundary  layer  are  also  obtained. 
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Returning  to  the  scaling  functions,  we  note  that  equation  19  requires 
do'/dR^  ,  and  hence  for  simplicity  we  have  treated  o'  as  one  of  the 
dependent  variables  of  the  final  system  of  ordinary  differential  equations. 
We  differentiate  equation  17  which  becomes 


.2  -2 


-■f-V  *  ‘  —  -  P-^r  <1  +  m  )  +  — f2l 

"  dR-  L0'2  “  ^  j  H»- 


dR- 

X 


+  [ <f>  [■x*-  -  l)  +  2<f2>  A  fl  +  2A<f2>  f2  -L 

L  V  iO00  /  J  A  D_  —  J— 

Combining  equations  (18)  and  (19) 


0 


dR- 
x 


u  dR— 
e  x 


(21) 


dRx  1 

and  recalling  ——  ■  — y  ,  it  follows  from  equation  (19)  that 

dR-  * 
x 


(22) 


—0  _  RT*  An  *  _  RT*  du 

f 2  -  H  -V  “  (1  +  H)  ~- 

”1  °  dR—  u  dR— 

x  ex 


dR- 
x 


f2  . 

(  T° .  J 

1  4.  _  u 

(R6* 

du 

e 

,  1  +  T  H 

i  W  1 

u 

e 

dR- 

X 

(23) 


Equations  (20) -(23)  and  the  integral  formulation  of  the  low  speed 
flow  (see  Appendix)  represent  the  total  formulation  which  is  composed  of 
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seven  ordinary  differential  equations  with  seven  dependent  variables; 

Rt*»  f,  u  /u^,  a’»  n  ’  ,  R  with  R—  as  the  independent  variable.  These 
equations  are  well  defined  provided  u  (R  )  is  specified,  and  the 

6  X 

variables  are  initialized. 

We  could,  of  course,  resort  to  the  particular  high  speed  experiment 
of  interest  for  such  initial  conditions.  However,  in  order  for  the 
procedure  to  be  useful  in  practice  (i.e.,  general  engineering  predictions) 
it  is  necessary  to  choose  these  values  independently  of  a  particular 
experiment.  We  have  chosen  the  values  tabulated  by  Coles ^  for  which 
the  wake  of  the  incompressible  boundary  layer  disappears  as  being 
characteristic  of  the  low  speed  turbulent  boundary  layer  near  its  "origin". 


Cf  -  5.9  x  10  J 
o 

R-  -  A. 25. 

U 

O 

At  the  "origin"  we  have  taken  R^  -  Rjf  -  0  and  (u  / u  )  *  1  ; 

Ao  o  e  °°  o 

a  '  can  be  calculated  from  equation  (17),  M„  =  M  and,  of  course, 

O  '  Cq  cxi  *  » 

n  '-o', 

o  o 

The  wall  temperature,  T  /Tn  ,  has  been  carried  throughout  the 
analysis  as  arbitrary  and  hence  formally  non-adiabatic  boundary  layers 
could  also  be  predicted.  In  a  strict  sense,  however,  Crocco's^^ 
suggestion  that  the  heat  transfer  of  the  two  flows  must  match  can  be 
shown  to  be  a  requirement  of  the  mapping  [see  Lewis  ].  In  particular 
we  must  have 


3T 


JL  (21) 

T  9y 


w 


and  hence  the  low  speed  flow,  while  incompressible,  cannot  be  of  constant 
density  if  it  is  to  correspond  to  a  non-adiabatic  compressible  flow. 
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4.  APPLICATION  OF  THEORY 


4.1  LOW  SPEED  -  CONSTANT  PRESSURE 

In  order  to  evaluate  the  formulation  for  the  low  speed  turbulent 

boundary  layer,  particularly  in  light  of  our  choice  of  initial  conditions, 

we  have  made  comparisons  with  the  experimental  data  of  Wieghardt^ 

—  __ 

(Figure  5).  C^  and  R^  are  seen  to  be  veil  predicted. 

4.2  SUPERSONIC  -  CONSTANT  PRESSURE 


Figure  6  shows  typical  comparisons  of  the  theory  and  experimental 
constant  pressure  adiabatic  data  of  Coles. ^  The  agreement  is  remarkably 
goo(  considering  that  no  recourse  was  made  to  the  high  speed  experiment 
for  the  initial  conditions  used  in  the  calculations. 

4.3  SUPERSONIC  -  VARIABLE  PRESSURE 


There  appears  to  be  a  remarkable  lack  of  consistent  detailed  experimentation 
on  two-dimensional,  supersonic  turbulent  boundary  layers  in  variable  pressure 
that  do  not  involve  separation. 

We  have  chosen  the  shock-wave  impingement  experiments  of  Hakkinen^1^ 
at  Mach  1.5  to  model  theoretically.  The  predicted  shear  stress  is  in 
agreement  with  the  experimental  values  (Figure  7)  and  the  characteristics 
of  the  equivalent  low  speed  flow  are  shown  in  Figure  8. 


In  order  to  explore  the  theoretical  predictions  at  higher  Mach  number, 
we  have  taken  a  gradually  increasing  pressure  C  ■  C  (R^)  at  various 
Mach  numbers.  The  results  shown  in  Figure  9  show  an  interesting  and 
somewhat  surprising  reversal  of  the  behavior  of  shear  stress  with  increasing 
Mach  number. 


^  1  — —  ^  C. 

H  was  found  to  be  in  good  agreement  with  the  experiment  beyond  Rx  ■  2  x  10  , 
but  there  exists  a  discrepancy  at  lower  Reynolds  numbers  which  is  believed 
due  to  the  neglect  of  the  laminar  sublayer  in  the  low  speed  integral 
formulation.  This  is  not  expected  to  be  an  important  limitation  for  most 
high  speed  turbulent  boundary  layer  situations. 
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This  result  as  well  as  the  theory  in  general  must,  of  course,  be 
verified  by  future  experimentation. 

4.4  BOUNDARY  LAYER  SEPARATION 

We  find  that  in  the  limit  as  f  ->  0  ,  the  low  speed  formulation  predicts 

that  (ti/u  )  (du  /dx)  ■  -  4.05  x  10  ^  . 
e  e 

By  recalling  equation  (12)  and  noting  that 


we  find  that,  independent  of  the  history  of  the  turbulent  boundary  layer, 


The  result  is  in  good  agreement  ^*10%)  with  the  empirical  correlation 
of  Zukoski^^  which  represents  a  compilation  of  adiabatic  wall,  turbulent 
boundary  layer  separation  data  for  -  2  ->  6  . 

4.5  SUPERSONIC  EQUILIBRIUM  PROFILES 

Clauser^^  has  shown  that  by  adjusting  the  pressure  distribution  such 

that 


<5  dp  — 

— - f  -  0  -  constant 

x  dx 
w 

the  outer  portion  of  low-speed  turbulent  boundary  layer  velocity  profiles 
can  be  maintained  in  a  state  of  "equilibrium",  i.e.,  they  are  invariant 
in  the  sense  that 
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By  following  the  rules  of  the  transformation  and  invoking  the  pressure 
gradient  relationship  (equation  12) ,  we  can  generalize  this  concept  to 
supersonic  flows.  The  supersonic  equivalent  of  Clauser's  low  speed 
equilibrium  boundary  layer  is  predicted  to  be  of  the  form 


where 


and 


(1  -  m  H) 
e 


* 

6  dp 

x  dx 
w 


The  modeling  outlined  in  the  paper  can  be  used  to  generate  the  required 
pressure  distributions  for  any  initial  boundary  layer  and  Mach  number 
corresponding  to  a  given  experimental  facility. 
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5.  CONCLUSIONS 


Summarizing,  we  have 

(1)  found  that  the  deviation  of  the  high  speed  velocitj  profiles 
from  the  equivalent  low  speed  profile  occurs  primarily  in  the  inner 
region  of  the  boundary  layer,  not  in  the  wake. 

(2)  given  evidence  that  the  discrepancy  is  due  to  increased 
"dissipation"  near  the  laminar  sublayer  in  the  compressible  flow  and 
concluded  that,  for  not  excessively  high  Mach  numbers,  this  should 
not  affect  the  application  of  the  theory  to  obtain  overall  results. 

(3)  derived  the  equations  which  govern  the  scaling  functions  for  a 
variable  pressure  in  terms  of  the  integral  properties  of  the  low  speed 
turbulent  boundary  layer. 

(4)  set  forth  a  complete  formulation  of  the  modeling  of  the  supersonic 
turbulent  boundary  layer  to  its  low  speed  equivalent  and  developed  a 
computer  program  to  implement  it. 

4 

(5)  applied  the  theory  to  a  variety  of  supersonic  turbulent  boundary 
layer  flows. 
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Figure  2.  Lilect  of  Arbitrarily  Increasing  Dissipation  in  I’ipe  Flow. 
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Figure  5.  Incompressible  Integral  Properties  -  Constant  Pressure. 
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Figure  6.  Supersonic  Integral  Properties  -  Constant  Pressure. 
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Figure  7.  Supersonic  Integral  Properties  -  Variable  Pressure. 


26 


0.878  0.911  9.939  0.972  1.004  1.014  1.026  1.040 

tfyXlO'6 

DISTANCE  FROM  ORIGIN 


Figure  8.  Low  Speed  Equivalent  of  Hakkinen* 
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APPENDIX 

LOW  SPEED  INTEGRAL  FORMULATION 


qq) 

The  formulation  [see  Alber  ]  used  consists  of  the  integrated  form 
of  the  momentum  equation 


H  +  Rj*  +  R^-*  ~  ~  -  f2  -  (2H  +  1)  — 


dR- 

x 


dR- 

x 


Bf  dR- 
x 


u  dR— 
e  x 


(24) 


and  its  first  moment 


—  ^RT*  zT  aT  aT  rT*  du 

+  -  3J-A--S. 

dR-  3tt  dR-  6  3f  dR-  D  u  dR- 

x  x  x  ex 


(25) 


and  the  differential  form  of  the  skin  friction  law 


dR^r*  — 

“3—  +  R7-*  [P]  ~  +  Rj*  [Q]  ~- 
dR—  dR—  dR— 

X  XX 


R7*  du 
6  e 


u  dR— 
e  x 


(26) 


where  tt  comes  from  Coles 
layer 


,(D 


analytical  form  of  the  turbulent  boundary 


y  u 

-3-1  +  C  +  £  W(y/6) 


k  ■  .40  (Kantian' s  constant) 


C  ■  5.10  and  u^  ■  V  x^/p  ,  and  H,  J,  [P],  [Q]  and  CD  are  specified 
functions  of  *tT  and  f  (see  Reference  18). 


1-1 


In  order  to  arrive  at  the  dependence  of  on  tt  and  f 

_  _  _  _  _  -- 

necessary  to  link  dT  -  (<$  /t  )  (dp/dx)  to 


and  here  we  have  modified  Alber's  formulation  to  agree  with  the 

(19) 

values  predicted  for  equilibrium  flows  by  Mellor  and  Gibson 
the  following  form  was  chosen 


ST  -  -  .852  -  .0596 


G  +  (G/5.90) 


,  it  is 


numerical 
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the  pressure  gradients  cf  the  two  flows  and  predicting  the  behavior  of 
the  integral  properties  of  the  low  speed  flow.  The  resulting  formula¬ 
tion  predicts  the  integral  properties  of  the  high  speed  turbulent 
boundary  layer  as  well  as  the  required  modeling  to  an  equivalent  low 
speed  flow.  Finally,  the  theory  is  applied  to  a  variety  of  situations. 
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